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Context. Narrow-band interference filters can be tuned toward shorter wavelengths by tilting them from the perpendicular to the 
optical axis. This can be used as a cheap alternative to real tunable filters, such as Fabry-Perot interferometers and Lyot filters. At the 
Swedish 1 -meter Solar Telescope, such a setup is used to scan through the blue wing of the Can H line. Because the filter is mounted 
in a converging beam, the incident angle varies over the pupil, which causes a variation of the transmission over the pupil, diff'erent for 
each wavelength within the passband. This causes broadening of the filter transmission profile and degradation of the image quality. 
Aims. We want to characterize the properties of our filter, at normal incidence as well as at diff'erent tilt angles. Knowing the broadened 
profile is important for the interpretation of the solar images. Compensating the images for the degrading eff'ects will improve the 
resolution and remove one source of image contrast degradation. In particular, we need to solve the latter problem for images that are 
also compensated for blurring caused by atmospheric turbulence. 

Methods. We simulate the process of image formation through a tilted interference filter in order to understand the eff'ects. We test 
the hypothesis that they are separable from the eff'ects of wavefront aberrations for the purpose of image deconvolution. We measure 
the filter transmission profile and the degrading PSF from calibration data. 

Results. We find that the filter transmission profile diff'ers significantly from the specifications. We demonstrate how to compensate for 
the image-degrading eff'ects. Because the filter tilt eff'ects indeed appear to be separable from wavefront aberrations in a useful way, 
this can be done in a final deconvolution, after standard image restoration with Multi-Frame Blind Deconvolution/Phase Diversity 
based methods. We illustrate the technique with real data. 
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1. Introduction 



Narrowband (NB) interference filters can be tuned in wavelength 
by varying the incident angle. This is often used to allow for 
filter manufacturing margins by specifying the peak to the red 
of the wanted wavelength and then mounting the filter at a small 
angle to compensate. At the Swedish 1-m Solar Telescope (SST; 
[Scharmer et al.||2003a ), the eflfect is used for scanning through 
the blue wing of the Can H line by use of a narrow-band (NB) 
filter mounted on a motor-controlled rotation stage. The tuning 
is possible because the filter cavity is geometrically wider for a 
ray intersecting with the filter at an angle from normal incidence. 

However, even when the filter is mounted with no tilt with 
respect to the optical axis, rays arrive at a variety of angles. This 
has consequences for the spec tral resolution as well as for the 
spatial resolution ( [Beckers 1199 8 ). Depending on the optical con- 
figuration, the spread in incident angles may cause a broadening 
of the transmission profile or variations in central wavelength 
over the field of view (FOV) or both. A problem related to the 
broadening is that the resolution suffers in telescopes that oper- 
ate near the diflTraction limit. The amplitude of the pupil trans- 
mission coefl&cient, as seen through the filter, is apodized by an 
amount that varies with the wavelength. Beckers evaluated the 
resolution eflTects of a Fabry-Perot interferometer (FPI) used in 
the pupil plane of a telecentric setup. | Von der Liihe & Kentischer 
( [2000] ) i ncluded also the phase of the field transmission in their 
analysis. [ Scharm er] ([2006 ) showed that the dominating part of 
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the phase eflTects is in the form of a wavelength independent, 
quadratic pupil phase that can be compensated with a simple re- 
focusing of the detector. 



The eflfects investigated by the authors cited above are 
greater the narrower the bandpass, making this a problem pri- 
marily for FPIs with FWHMs in the picometer range. However, 
as we will demonstrate in this paper, there are significant eflfects 
also for filters with a FWHM of ~0. 1 nm when tilted by a few 
degrees. We analyze a NB double-cavity filter, mounted in a con- 
verging beam at an angle from the optical axis. We show that the 
resulting degradation of image quality can be compensated for 
by deconvolution and that images degraded by this effect as well 
as atmospheric seeing can be significantly improved. The effect 
is apparently separable from phase aberrations in the sense that 
they can be modeled as separate transfer functions, to be applied 
in sequence. This is beneficial for wavefront sensing with the 
Multi-Frame Blind Deconvolution/ Joint Phase Diverse Speckle 
(MFBD/JPDS) formalism of [Lofdahl| ( [20021 ). 



We illustrate the technique with images from the SST 
recorded through a Can H filter mounted on a computer con- 
trolled rotating stage. We use simulations of such data to illus- 
trate the separability of the apodisation effects and phase aber- 
rations. Appropriate and sufficient information needed to com- 
pensate for the image degradation can be obtained by calibration 
using pinhole images. 
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2. Theory and notation 

2.1. Tilted cavities 

In order to understand what happens when we tilt a filter, we be- 
gin by deriving the transmission coefi&cient for the electric field. 
The matrix formalism described byl Klein & Fur tak (1986, Sect. 
5.4) implements multiple reflections and transmissions in both 
directions through many diflferent layers. Assume light incident 
on a stack of N layers from the left. We define the column vec- 
tors 



and 



(1) 



where Ej represents the electric field on the left side of layer j 
and E^. the field on the right side of the same layer. Indices L and 
R denote field moving to the left and right, respectively. 

The relationship between the field in the incident medium 
and the field in the final medium can then be written as 



E^ - Sia^Ea^, 

with the stack matrix 



^11 ^12 
^21 ^22 



H12L2 • • -La^-iHa 



(2) 



(3) 



The interface transition matrix Hij and the layer matrix are 
defined as 



1 Pij 
Pij 1 



and 



e^^ 



(4) 



repectively, where Ttj and ptj are the transmission and reflection 
coeflftcients, respectively, of the interface between medium / and 
medium j, and 



47r , 
^7 = y^A-cosay. 



(5) 



Here rij is the refractive index of medium j, hj is the geometrical 
thickness of the cavity, and aj is the angle of the ray from the in- 
terface surface normal within medium j. Using these definitions 
and the boundary condition Ei^^ = (in the final medium there 
is no field moving to the left), the transmission coefl&cient of the 
stack can be written as 



^Rl 



1 

^22 



(6) 



Applying this formalism to a single cavity with index ^2 be- 
tween media of index rii = ^a^, we get 



Tl2T21^' 



-^2 



If n2, h2, d2, T12, T21, P12, and pu are known, like for a FPI, this 
expression can be used to calculate the transmission of the cav- 
ity for any incident angle. However, manufacturers of dielectric 
filters usually only disclose the intensity transmittance (T = |tP) 
at normal incidence. 

The transmission coefliicient is a function of the wave- 
length as well as the incident angle through ^2- When the an- 
gle changes, the wavelength needed to keep 132 (and thereby r) 
constant also changes. The effect of changing the angle is there- 
fore to shift the transmission profile in wavelength. With the as- 
sumption of small angles, repeated application of Snell's law. 



and trivial algebra, the well known (e.g., Smith|1990 chapter 7) 
expression 



=^^/l-l^sinal , 

\n2 



(8) 



can be derived, where is the shifted wavelength and we define 
a = ai. 

With double cavities we get 



T12T23T34 • e 



1 -\-p^e-'^^ -p23p34e-'^^' -Pi2p23e-'^^' ' 



(9) 



where we have defined (5 = {f52+ and p = P12P23 + P23P34 + 
P12P34. The numerator shifts as given by the shift of which 
means it behaves like the transmission profile for a single cavity 
and an "eflfective" refractive index n can be defined. The last two 
terms in the denominator shift diff"erently but the behavior of the 
numerator should dominate. 

2.2. Pupil transmission in a converging beam 

For a piece of optics in a converging beam, the incident angle 
of a light ray varies with the orientation of the optics relative 
to the optical axis, but also with the location within the pupil 
that the ray emanates from. As we have seen, the transmission 
coeflftcient t(A) of an interference filter varies with the incident 
angle, so that it is shifted toward shorter wavelengths. 

For incoming light with a wavelength A, the generalized 
pupil function may be written as 

P(x^,y^,A;(p,e) = Tp(Xp,};p,/l;^)A(xp,};p)exp{i0(xp,3;p)}, (10) 

where (Xp, jp) are coordinates within the pupil (here normalized 
to unit diameter), A is a binary function representing the geo- 
metrical shape of the pupil stop, i is the imaginary unit, 6 is the 
tilt angle of the filter around the y axis, and is the phase aber- 
rations from the atmosphere, the telescope, and any other optics, 
independent of A. The filter's complex transmission coeflftcient 
as seen through the pupil can be written as 



Tp(xp,};p,/l;6>) = Tf{A\A,a(x^,y^;0))), 



(11) 



where Tf is the normal incidence complex transmission coeflft- 
cient of the electric field vector through the filter, and A' is the 
shifted wavelength given by Eq. ([8]). In Appendix [A| we derive 
an expression for the incident angle, a, on a tilted filter in a con- 
verging beam. With the added assumption that the detector is 
small, we can write it as 



(7) Q'(xp,3;p; ^) = arccos 



F cos 6 -\- Xp sin 



(12) 



JF2 + X^+j2)l/2j' 

where F is the effective F-ratio. 

2.3. Point spread function and optical transfer function 

Because of Eq. ([8|, we know that NB interference filters can be 
tuned by tilting the filter normal from the optical axis. In ad- 
dition to the desired change in peak wavelength, the tilting has 
consequences for the width of the profile and for the shape of the 
PSF because of Eq. ([12]). 

We model the image formation process as a space-invariant 
(valid for sufl&ciently small sub-fields) convolution between an 



2 



Mats G. Lofdahl et al.: A tilted interference filter in a converging beam 



unknown object / and a point spread function Sk, the Fourier 
transforms of which are F and S k, where the index k corresponds 
to a particular data frame. An additive Gaussian noise term rik is 
assumed. The Fourier transform Dk of the observed image dk is 
then related to F, Sk and Nk via 



360 



Dk = FSk + Nk 



(13) 



These quantities are all functions of the spatial frequency coor- 
dinates u and V but to ease the notation, we will not write them 
out. The assumption of additive Gaussian noise leads to a simple 
Maximum Likelihood expression for an estimate of the object 
through multi-frame deconvolution, 



F = HY,DkSllY,\S 



(14) 



where H is a low-pass noise filter and ^ used as a superscript 
denotes the complex conjugate. 

Because the profiles are shifted by diflferent amounts over 
the pupil, the transmittance of light of different wavelengths 
through the pupil is not uniform. At every wavelength within 
the broadened passband, the transmittance of light from diflferent 
parts of the pupil corresponds to different parts of the normal- 
incidence filter profile. Because light interferes only with light 
of the same wavelength, S k cannot be calculated with the usual 
quasi-monochromatic assumption. Instead, it must be modeled 
as the integrated, monochromatic contributions from the differ- 
ent wavelengths within the passband. With a discrete approxi- 
mation, we can write 



Sk = S((Pk, e) = SAY,^ {|r ' P(x^.yp, h\ 0^, ^)| 



(15) 



where the summation is performed over L samples distributed 
over the broadened passband. 

What we would really want is for 5" to be separable in two 
parts: 1) the usual OTF, including the pupil geometry, the atmo- 
spheric wavefront, as well as any other wavefront components 
introduced in the optics, and 2) another transfer function, com- 
pensating for the filter tilt effects. We would like to write 



where 

^wf(0) = A(xp,jp)exp{i0(xp,jp)j 
and 

1 ^ 

^ tilt(^) = 7^^^ / - ^ mono(^^, ^), 

where we define the monochromatic tilt OTF as 



(16) 



(17) 



(18) 



(19) 



In the expression for 5'tiit, we choose to include the pupil ge- 
ometry and then divide it away as the diffraction limited MTF, 
^0 = ^wf(O)- This way it becomes natural not to define Tp out- 
side of the pupil. For simplicity of notation, and for the purpose 
of this discussion, we incorporate also any intentional phase di- 
versity (PD) focus contribution in 0^^. 




We explore this separability in Sect. |3. 3 [ below. The advan- 
tages with separability, for our purposes, is that the correction 



396.6 396.8 397.0 397.2 

Fig. 1: Filter profiles for normal incidence on the Can H filter. 
Solid black: Transmittance, Ff^X), measured by the manufacturer 
before delivery of the filter; Dashed red: theoretical transmit- 
tance; Blue: theoretical transmission phase, arg(Tf(/l)). 



can be applied after restoration of the images for atmospheric 
turbulence effects, i.e., MFBD-based methods or speckle inter- 
ferometry. And, since we assume the variation over the field of 
view is negligible, it can be done for the entire FOV in a sin- 
gle operation instead of by subfields like for the anisoplanatic 
atmosphere. 



3. A 0.1 nm wide Ca ii H 396.9 nm filter 

3A. Filter properties 

As the specific case, motivating this research, we discuss a 
Can H filter, regularly used at the SST. The filter, manufactured 
by Barr Associates, Inc., was delivered in August 2000. The 
coating parameters are not public so we cannot calculate Tf but 
along with the filter we got the measurements shown in Fig. T] 
The central wavelength was specified to Ac = 396.88 nm, the 
peak transmittance, max Ff, to 44%, and the FWHM to 0.1 1 nm 
(at 23 °C and normal incidence). 

In order to use Eq. ([8]), we need the effective refractive index. 
Rouppe van der Voort| ( 20021. page 26) contacted the manufac- 
turer after the filter was delivered and got the value n = 1.596. 
Because we wanted to explore the effects of the phase, |Potter| 
( 2004) later provided some design data: the phase of the filter 
transmission coefficient, arg{Tf(/l)}, as well as the curve form of 
the transmittance profile in digital form, see the red and blue 
curves in Fig. [T] The wavelength range for both quantities is 
394.9 nm < A < 398.9 nm. We show here the central part, 
rescaled to the measured peak transmittance of 44%. 

The light beam from t he telescope forms a pupil image, 
where our adaptive optics ( Scharmer et al.| |2003b ) deformable 
mirror is mounted. From there, the beam is re-imaged by a field 
lens with a focal length of ~ 1.5 m, making a F/46 beam. As seen 
from the pupil image, through the re-imaging lens, the 7 mm 
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Fig. 2: The transmission coefficient, Tp, profile data cube, cal- 
culated from t he m anufacturer's theoretical profiles by use of 
Eqs. (p^-([r2|).|(a)]The modulus f or = 0, ItJx^^O, A;0)l su- 
perimposecHDr & = 0°, 3°, and 6° as labeled. [(b)| The pupil apodi- 
sations, |Tp(Xp, jp, /l^; 6°)|, sampled at the wavelengths indi- 
cated with the orange lines in |(a) (c) The corresponding pupil 
phases, arg{Tp(Xp, jp, A^; 6°)}. 



center to comer distance of our 10 mm by 10 mm detector cor- 
responds to an angular size of 0.007/1.5 rad or (f.21. As we will 
see, angles of this size has negligible impact so we are justified in 
ignoring variations over the FOV. The image scale, O'.'034/pixel, 
made by the re-imaging le ns was establish ed by observations of 
the Venus transit in 2004 ( Kiselman|2008| ). 

The pupil apodisations vary with tilt angle, which is illus- 
trated in Fig. [2] Increasing the angle moves the passband from 
the Ca II H line core (formed in the solar chromosphere) toward 
the blue wing (formed in the photosphere). In the axis corre- 
sponding to the tilt direction, the transmittance profile peak po- 
sition varies over the pupil (Fig. [2a]). In the perpendicular direc- 
tion, the peak position varies only slightly with distance from the 
center of the pupil. Sample pupil apodisations for one of the tilt 
angles are shown in Fi g [25] Because of the angles of the trans- 
mittance ridges in Fig. |2a[ larger tilt angles result in narrower 
apodisations over the pupil for each wavelength. For 0?0, the 
transmittance varies with wavelength, but with small variations 
within the pupil. For larger tilt angles, the apodisation function 
looks more and more like a band sweeping across the pupil. 

The transmittance phase, arg(Tf), has an approximately lin- 
ear gradient where the modulus, |Tf|, has its peak, see Fig.[T] The 
phase profiles are shifted the same way as the transmittance pro- 
files, so the phase gradient peaks at the same location in the pupil 



as the transmittance, see Figs.pmand 2c 



3.2. Profiles 

Integrating over the pupil results in broadened profiles, as noted 
by Beckers (1998]), because of the varying shifts in the Xp di- 
rection. For our lilter, ^ = 0?0 results in a transmission profile 
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Fig. 3: Integrated filter transmittance profiles, Zx^^j^ ^p(^' ^) 
(red) corresponding to filter tilt angles (6) as indicated in the 
figure. The ^ = 0° profile is virtually identical to the normal in- 
cidence filter profile, Tf(A) (blue). The black line is the Can H 
spectrum of the synthetic data of Sect. |4] 



that is virtually identical to the normal incidence profile but for 
angles of a few degrees, the broadening is significant, see the 
broadening as calculated for the manufacturer's profile in Fig. [3] 
In order to test the predicted broadening and wavelength 
shifts, the filter was examined using the TRI-Port Polarimetric 
Echelle-Littrow spectrograph (Kiselman et al. in prep.) at the 
SST on 2008-05-13. The filter on its rotatable mount was placed 
in front of the spectrograph slit and the telescope was pointed to 
solar disk center. Spectra (2000 wavelength points in the range 
396.1 nm < A < 397.1 nm) were collected without the filter, as 
well as with the filter at 18 diflferent tilt angles that turned out 
to be -2? 16 < < 5?84, i.e., angles on both sides of the sym- 
metry point are well represented. The spectra were reduced with 
standard methods. Filter profiles were then calculated by nor- 
malizing the filtered spectra with the unfiltered spectrum. The 
resulting profiles contain small residues from blending spectral 
lines which have not been completely removed by the normaliza- 
tion. This is due to the presence of stray light within the spectro- 
graph, which is corrected for in the reductions but the procedure 
obviously does not work perfectly for spectra of this kind. We 
deem the resulting profiles good enough for our current purpose, 
however. 

The wavelength scale was calibrated with the help of a spec- 
tral atlas (Bra ult & Neckel|1987j and transformed to the labora- 
tory rest frame. The zero angle is then given by symmetry and 
the peak wavelength for normal incidence, A^ = 396.89 nm, fol- 
lows. 

In Fig. [4] we show the transmittance profile measured with 
the smallest angle, 0?16. For such small angles, while tilting 
shifts the central wavelength, it does not broaden the profile sig- 
nificantly. Therefore this should be a good measure of the normal 
incidence profile Tf. We demonstrate that this is so by relocating 
the measured profile to the symmetry point. Using this to model 
the measured curve with the filter tilt and pupil apodisation for- 
malism, it is clear that this is just the amount needed for the tilt to 
shift it back to the measured position and that the match in pro- 
file shape is excellent. The only thing that happens to the profile 
shape is that the small artifact in the red flank is smoothed out 
slightly. The FWHM of the measured Tf is 0.097 nm rather than 
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TRIPPEL, measured @ 0=O?16, FWHM=0.097 nm 
TRIPPEL, "normal incidence" 
TRIPPEL, "normal incidence", shifted to 0=O?16 
Theoretical, normal incidence 
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Fig. 4: The filter normal incidence transmittance as measured 
with TRIPPEL. Green: the measured profile shifted towards the 
red to compensate for tilt angle, our measured Tf. Red: model 
of black curve using measured Tf. Black: measured profile with 
smallest angle. Blue dashed: theoretical Tf (peak reduced to 
42%). 



TRIPPEL, "normal incidence" 
TRIPPEL, measured 
TRIPPEL, shifted 
Theoretical, normal 



and shifted 




> i_5?84 



396.2 396.4 396.6 396.8 397.0 

X / ^ nm 

Fig. 5: TRIPPEL transmittance measurements at ^ > 2°. Black: 
measured profiles. Red: models of measured curves using the 
TRIPPEL-measured "normal incidence" profile. Blue dashed: 
models of measured curves using the theoretical normal inci- 
dence profile with peak reduced to 42%. 



the 0.11 nm value from the manufacturer, a significant diflTerence 
of more than 10%. 

The consequences for large angles are shown in Fig. [5] Using 
the TRIPPEL-measured normal-incidence profile, we can model 
the measurements very well, which means the profile shift and 
broadening based on single-cavity theory work well and that the 
eflTective refractive index, n = 1.596, must be a good approxima- 
tion. Shifting the manufacturer's profile produces profiles that 
broaden in the right way, but where the peak transmittance does 
not decrease with tilt angle fast enough. This is caused by the 
wider normal incidence profile and the fact that the broadening 
preserves the area under the profiles. 

If a manufacturer is not able or willing to provide an eflTective 
refractive index, it should be possible to infer a useful estimate 
from data like these, by fitting Eq. ^ to the peaks of the shifted 
profiles. 
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Fig. 6: Top: Diflfraction limited MTFs modified by the filter tilt 
eflTects, • l^tiit(^)l; Bottom: Tilt MTFs, |^tiit(^)|. Tilt angles 6 
as indicated. 




(a)OPO (b)l?5 (c)3?0 (d) 4?5 (e) 6?0 (f)T5 

Fig. 7: Top: PSFs corresponding to • ^tiit(6>); Bottom: PSFs 
calculated with arg{Tp} set to zero. Tilt angles 6 as indicated. 



3.3. Point spread functions and Strehl ratios 

The pupil phase for any particular varies with A, so it is not 
obvious what phase, if any, would be eflfective as a correction of 
the summed 5'tiit. In this section we investigate this. We use the 
theoretical Tf profile here, because we have modulus and phase 
that go together. 

The computed MTFs shown in Fig. [6] demonstrate that the 
tilt is a one-dimensional effect; there is no eflTect at all in the y 
direction, perpendicular to the tilt. 

The PSFs corresponding to large angles shift significantly in 
the X direction, as can be seen in Fig. [7] This indicates that there 
is a significant common tilt in the x direction in the monochro- 
matic wavefronts. This can be understood by noting that in 
Fig.[2j the monochromatic phases have gradients in the same po- 
sitions that there are peaks in the modulus. Ignoring the phases 
produces PSFs that are not shifted and also do not contain the 
asymmetry in shape that is also evident in the figure. 

In real data, we see much larger shifts than the few pixels 
visible in Fig. [7] In fact, most of the shift is caused by refraction 
through the glass substrate. An order of magnitude estimate of 
this shift can be calculated by use of the formula d = t6(n -I) In 
given by| Smith| p990, Chapter 4). With n = 1.596, thickness 
t ^ 5 mm, and tilt angle = 6?0 = 0.10 rad, we get a shift of 
d ^ 0.20 mm = 26 pixels of size 7.4 jim. 

Scharmer] ( |2006 ) found an optimum focus compensation for 



the phase eflTect of a non-tilted FPI by maximizing the Strehl 
ratio. In Fig. [8] we show Strehl ratios calculated with and with- 
out taking the phase of the transmission into account. The Strehl 
ratio without phase is reduced from almost unity for ^ = to 
-0.6 for 7? 5. Including the phase reduces the Strehl further, par- 
ticularly for larger 6. Maximizing the Strehl ratio by varying a 
common focus compensation we find negligible improvements 
in Strehl ratio. A tenth of a percent for the small angles, even 
less for the larger angles. The wavefront RMS of the estimated 
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e/1 

Fig. 8: Theoretical Strehl ratios for the Can H filter, calculated 
by integrating the MTFs. The pupil functions at each wavelength 
are calculated without the phase, with the phase, and with two 
diff'erent compensations for the phase. O4 is an optimum fo- 
cus term, while On is an optimized sum of focus, astigmatism, 
coma, trifoil, and spherical aberrations. 



focus terms ranges from 0.010 waves for 6 = 0?0 to 0.006 waves 
for6> = 7?5. 

The individual apodisations and phases for the non-zero an- 
gles are not circularly symmetrical so astigmatism and coma 
would not be unexpected. However, including Zernike poly- 
nomials 4-11 in the optimization gives insignificant additional 
Strehl improvements and no change in the estimated wavefront 
RMS [^Changes, if any, are in the sixth digit. There is no reason 
to expect higher order modes to have a significant eff'ect. 

The compensated Strehl ratios are also plotted in Fig. [8] 
where Om = Yjm=4^m^m are the optimized phases. are 
Zernike polynomials, m = 4 corresponds to focus. It is apparent 
that the compensations do not alter the Strehl ratios significantly. 

This suggests that, at least for this particular setup, the ef- 
fect of the varying incident angles on the filter should inter- 
fere very little with our normal Multi- Object Multi-Frame Blind 
Deconvolution (MOMFBD; [van Noort et al.||2005| ) restoration 
practice. Image shifts are already routinely taken care of in cali- 
bration of camera alignment by use of a pinhole array. The focus 
term is compensated already when the camera is focused, usu- 
ally for 6 = 0?0. This leaves an uncompensated focus of only 
0.004 waves or less at the larger angles, which is insignificant 
compared to the accuracy of our camera focusing. Evident from 
the way the PSF moves with angle when the phase is included 
is that there is also a global tilt term that varies with 6, but this 
does not change the shape of the PSF. 

This means the OTF is to a good approximation separable 
into an atmospheric wavefront part and a correction for the filter 
tilt as in Eqs. ([T6])-([T8]). The smearing caused by the elongated 
PSFs aff'ects the SNR by lowering the MTF but should otherwise 
pass through the MOMFBD processing without changing the es- 
timated phases. The restored object is a version of the real object 
that is convolved with the PSF corresponding to the apodisation 
eff'ects without the wavefront phase. This can be taken care of 
with a simple post-restoration deconvoluti on st ep. The conse- 
quences for alignment are discussed in Sect. |5. 3 [below. 



4. Synthetic solar data and expected errors 

The Strehl experiments correspond to objects that are point-like 
and do not depend on A. However, the solar structures vary with 
wavelength within the passband of the filter, particularly if it's 
broadened and particularly where the line gradient is large. In 
this section we simulate the formation of solar images by sum- 
ming the contributions from many wavelengths in a synthetic 
data cube. We make "perfect" images as well as images based on 
convolution with either 5*0 or 5'mono before summing. We evalu- 
ate how well the "perfect" data can be recreated by deconvolving 
the summed images. 

The synthetic image data cube was calculated from 3D MHD 
simulati ons kindly provided by Ma ts Carlsson and bas ed on the 
code by [Stein & Nordlund[ 09^ . We used MULTI ( [Carlsson[ 
1986[ ), a Can H 5-level atom model, and line opacities from 
the V ienna Atomic Line Database ( [Piskunov et al .' 1995 ; Kupka| 
jet al. 1999 , 2000) for the line blends. The data cube comprises a 
wavelength range covering the tunable filter range (394.2 nm < 
A < 399.5 nm, SA = 0.01 nm). The most relevant part of the 
spectrum is shown in Fig. [3] 

Because we have established separability in Sect. [3.3| we do 
not have to include atmospheric turbulence in these simulations. 
The experiment corresponds to images that are perfectly restored 
from those eff'ects by image restoration with methods such as 
MOMFBD or Speckle interferometry. We included neither noise 
nor stray light in these simulations. 

The results of the experiments can be found in Fig. [9] The 
RMS errors are calculated as the standard deviation of the diff'er- 
ence of an image and the "perfect" image, after subpixel align- 
ment. The errors therefore measure only the residual blurring 
and not errors from image shifts caused by off'-center and asym- 
metrical PSFs. 



^ We used Zernike polynomials ordered as specified by'Nol?(fl976]), 
where indices 4-1 1 correspond to focus, astigmatism, coma, trifoil, and 
spherical aberrations. For the ID optimization of the focus term alone, 
we used Brent's Method. Starting from the focus-only optimum, we 
used the Downhill Simplex M ethod (AMOEBA) for Zernikes 4-11. For 
both methods, see |Press et al.| ( |1986j ). 



The circles in Fig. 9a represent the RMS contrasts of the 
"true" images. These contrasts vary between 23.5% and 33.3% 
because of the varying heights in the solar atmosphere where 
they are formed. Convolving each monochromatic image in the 
synthetic data cube with the diff'raction-limited PSF before sum- 
ming, i.e., multiplying the Fourier transform with 5*0, makes the 
contrasts drop as shown with the blue crosses, resulting in RMS 
errors in the range 4-5% as shown in Fig. [9b] Deconvolving with 
5*0 recreates the original contrasts and brings the errors down to 
zero as expected (blue pluses). 

For the remaining experiments, 5'niono was applied to each 
monochromatic image before summing. The contrast drops a bit 
more than for 5*0, particularly for the larger tilt angles and the 
RMS errors are larger (green crosses). Consequently, deconvolv- 
ing with 5*0 works well for the small angles but not for the larger 
angles (green pluses). In fact, for (f.O and 1?5, it does not matter 
much whether we deconvolve with 5*0, 5'tiit or any of the imper- 
fect OTFs described in the next paragraph. 

For the larger angles, deconvolving with 5'o-5'tiit also restores 
the original contrasts and reduces the errors to small fractions 
of the S 

mono versions. However, when using imperfect versions 
of 5'tiit (red symbols) the results for the larger angles is worse. 
Not using the phase of the transmittance profile, arg(Tf), gives 
approximately the same RMS errors as correcting only for the 
MTF, 5*0 • 15'tiitl, which is expected since neither can represent 
the asymmetry of the PSFs (red upward and downward triangles, 
resp.). Even worse is to use the asymmetric PSFs, but reversed 
(red diamonds). 
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Fig. 9: RMS contrasts (a) and RMS errors (b) from synthetic image deconvolution experiment, both in percent of the mean intensity. 
Black circles (O, in a only): Synthetic "true" images. Blue symbols: monochromatic images convolved with 5*0 and then summed. 
The crosses (x) correspond to those images and the pluses (+) to the same images after deconvolution with Sq. Green symbols: 
monochromatic images convolved with 5' mono (^.^)- The crosses (x) correspond to those images and the pluses (+) to the same 
images after deconvolution with Sq. The squares (□) correspond to those summed images again, but now deconvolved with • 
StiitiO). Red symbols: convolved images like for the green symbols but with three types of imperfections in the deconvolution. Up 
triangles (A): transmission phase not used for calculating 5'tiit. Down triangles (V): deconvolve with the MTF, Sq- |5'tiit(^)|. Diamonds 
(O): deconvolution with • 5'tiit(^) mirrored in the x direction. 



5. Real data and consequences for image 
restoration 

In this section, we work with SST data collected by Henriques 
et al. (in prep.) in May 2010. We used four MegaPlus II es4020 
cameras with 2048x2048 7.4-jim pixels. The beam was F/46 and 
the image scale is O'.'034/pixel. There was 23% oversampling in 
the Fourier domain. 

One camera collected data throught the tilt-mounted NB fil- 
ter. The observing programme specified tilt angles = 1?2, 2?5, 
3?0, 3?7, 4?2, 4?8, and 6?4. These angles correspond to pass- 
bands with the central wavelengths Ac = 396.84, 396.74, 396.67, 
396.57, 396.47, 396.34, and 395.93 nm, respectively, chosen to 
sample interesting heights in the solar atmosphere. In addition, 
we had two cameras behind a wideband (WB) filter centered on 
395.37 nm (FWHM 1.0 nm), one in the conventional focus and 
one defocused by 7 mm for approximately 1 wave peak-to-peak 
of focus PD. The fourth camera was mounted behind a fixed NB 
(FWHM 0.1 nm) 396.47 nm wing filter. 

We routinely perform a calibration step, where we collect 
data while scanning the filter through a large range of angles, 
approximately centered on ^ = 0° (in the case of these data: 20° 
with 0?1 resolution). We calculate the average intensities of the 
images and find the zero position by looking for the symmetry 
point. Because these data depend on the unknown spectrunj^ we 
do not model these data in order to draw any conclusions besides 
the zero tilt angle. 

5.1. Pinhole images and OTF correction 

As part of the standard calibration procedures at the SST, we 
collect data through artificial targets mounted in the primary 
(Schupmann) focus. One such target, designed for measuring 
straylight, consists of six holes of diameters that range from 
20 jim to 1 mm. We can use these to test our calculated 5'tiit cor- 



^ The Liege and FTS atlases differ and neither is guaranteed to match 
the data at our position. 



rections. We use a 5 12x5 12-pixel subfield with only the small- 
est, barely resolved pinhole. 

In order to gain SNR, we co-added many frames (each with 
an exposure time of 10 ms), aligned to subpixel precision. The 
available number of frames from this campaign is 500 in each 
of the tilted NB wavelengths and 3500 in the WB. The indi- 
vidual frames are dark corrected by subtraction of an average 
dark frame. After co-adding, an extra dark correction step is per- 
formed. All summed images were normalized to the average in- 
tensity within the 1-mm hole. A refinement dark level for the WB 
image was calculated as the median of the outer few rows and 
columns of the 5 12x5 12-pixel subfield and subtracted from the 
WB image. The NB images were then dark corrected by subtrac- 
tion of the amount needed to make the total intensity the same in 
all images. 

The power of the pinhole image is isotropic in the WB but 
not in the NB because of the tilt effect. In the y direction the WB 
and NB powers are equal but in the x direction, the NB power is 
attenuated. We demonstrate this for 6 = 6?4 in Fig.[TO| compare 
the black and red curves. Note the dip in power at approximately 
50% of the diflfraction limit. There is a corresponding dip in the 
NB power, although not as clearly visible because the power 
spectrum is so noisy at the spatial frequencies beyond the dip. 
These dips are caused by a zero-crossing in the Fourier trans- 
form of the target. The Fourier transform of a pillbox function 
(representing a circular pinhole) is a radial jinq^ function. Such 
zero crossings appear closer to the origin the larger the hole is, 
so the smaller the hole, the better. In the limit of zero diameter 
we have a S function, the transform of which is a constant. 

Correcting the NB power by use of the theoretical 5'tiit is 
not enough, as evidenced by the blue curve. We show here 
the theoretical correction based on the profile measured with 
TRIPPEL. Using the manufacturer's theoretical profile results 
in even slightly less correction. 



^ Similar to a sine function, jinc(r) - /i(r)/r, where Ji is a Bessel 
function of the first kind, oscillates with decreasing amplitude but with- 
out a defined period. 
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Fig. 10: Stray light target pinhole power spectra in the x direc- 
tion, affected by the 6?4 tilt. Angular averages within ±22.5° 
wide sectors centered on the axis directions. Calculated using 
a 1024X 1024-pixel subfield of the detector FOV. NB: original, 
corrected with theoretical MTF, and corrected with measured 
MTF (or OTF, the corrected powers are the same). WB for com- 
parison (covered by the NB compensated line). The TRIPPEL 
measured profile was used for the theoretical correction. Using 
the theoretical profile gives even less correction. 




linn 



(a) 6?4 



(b) 3?0 



(c) 2°5 



(d) r.2 



Fig. 11: Measured MTF correction |5'tiit(^)| based on straylight 
target pinhole power. Top: Dirty; Bottom: Cleaned. Tilt angles 6 
as indicated. The circumferences represent the diflfraction limit. 
The orientation of the ridge is ~2° from the y axis. Compare with 
Fig.|6l 

However, because the pinhole is the same in WB and NB, it 
should be possible to measure 5'tiit. We can get a first estimate 
by taking the ratio of the Fourier transforms of the two images. 



tilt 



FS 



►tilt, 



(20) 



wf 



where 5'wf is here the OTF representing phase aberrations in the 
optics on the optical table. The results for a few tilt angles are 
shown in the top row of Fig.[TT] These "dirty" measurements are 
dominated by noise at the higher spatial frequencies. However, 
since we know from Fig. [6] that they should be constant in the y 
direction, it should be possible to clean them. 

The ring at about 50% of the diffraction limit corresponds to 
the jinc dip in WB power. We mask this ring-shaped zero cross- 
ing artifact and calculate the median along the direction that is 
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Fig. 12: MTF corrections, |5'tiit(^)|, in the x direction for 6 = 
r.2, 3?0, 4?8, 6?4 (from top to bottom). Blue: theoretical MTFs; 
Red: measured MTFs. 



supposed to be constant. We smooth the result and then construct 
a cleaned 2D 5'tiit. In the real data, this direction is not exactly 
parallel to the y direction, so we first find the orientation of the 
ridge, rotate 5'tiit, do the cleaning and rotate the result back to the 



original orientation. The result is in the bottom row of Fig. 1 1 



We show the theoretical and measured corrections for a few 



angles in Fig. 12 Note the discrepancies between the blue curves 
and the corresponding red curves ! For the present data, the limit 
to where we can measure 5'tiit in the x direction is set by the SNR 
in the WB pinhole image in the denominator to a little less than 
80% of the diflfraction limit. We define 5'tiit to be unity outside 
this limit. Noise in the NB (numerator) is worse for the smaller 
angles because the Sun is darker in the core of the spectral line, 
the irregularities caused by this are apparent in the figure. 

It is apparent that we need more exposures to construct 5tiit 
functions that go all the way to the diflfraction limit in the x direc- 
tion, and that are not noise dominated above 50% of the diflTrac- 
tion limit. The noisy NB bump in Fig. [T0| suggests that we need 
to increase the number of NB frames by an order of magnitude 
or more. With an image rate of 10 frames/s, it should be possible 
to collect -5000 images in 8 min. With seven tilt angles, this cor- 
responds to almost an hour. However, this kind of data would not 
have to be collected each day or even by each observer interested 
in this kind of correction. 



5.2. Solar images 

On 2010-05-23, we collected solar images with an exposure time 
of 10 ms near disc center (ji ^ cos 20° ^ 0.94). Together with 
the WB, WB PD, and fixed NB data, the tilt-tuned NB images 
were restored for atmospheric turbulence eflTects with MOMFBD 
using the 36 most significant atmospheric Karhunen-Loeve 
modes. Observations and MOMFBD processing will be de- 
scribed in detail by Henriques et al. (in prep.). 



Figure 13 shows power spectra for 6 = 6° 4. For this an- 
gle, the central wavelength of the NB passband has shifted to 
395.93 nm, well within the passband of the WB filter. Here, the 
granulation has the same power in the y direction in NB as in 
WB, see Fig. |13b| We can expect this to be the case also in the x 
direction, which can be used for testing our 5'tiit compensation. 
In Fig. 13a the MOMFBD restored NB power (red) is attenu- 
ated compared to the WB power (black). Just as for the pinholes, 
the theoretical 5'tiit compensation (blue) does not fully correct 
the asymmetry in power, but the measured compensation does 
(green). 

There are visible eflTects in the images, that correspond to the 
correction of the anisotropy in the power spectrum. Compare the 
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Fig. 13: Solar power spe ctra, angular averages within ±22.5° wide sectors centered on the axis directions, (a) In the x direction, 
where the tilt effects are; I (b) I In the y direction, unaffected by the tilt eflTects. The dashed line represents noise filtered data. 




(c) 

Fig. 14: An IT^xlT' solar image collected with 6?4 filter tilt, (a) Restored with MOMFBD; (b) Compensated for tilt efiTects by use 
of the measured 5'tiit; (c) Compensated and noise filtered. 





MOMFBD restored image in Fig. 14a and the same image after 
correction for 5'tiit in Fig. |14b[ Because 5'tiit makes point-like ob- 
jects elongated in the x direction, the eflfect is most easily seen 
in the smallest structures. For example, strings of small bright 
features oriented in the x direction are more clearly separated 
after 5'tiit correction, while such strings oriented in the perpen- 
dicular direction mainly get higher contrast with respect to the 
surrounding dark lanes. 

The 5'tiit compensation enhances the noise somewhat, par- 
ticularly for the larger angles, so we need to use a noise fil- 
ter. Figure |14c| shows the compensated image filtered with the 
noise filter in Fig. |15a[ The noise filter is very asymmetrical, as 
it should because the attenuated power in the x direction will fall 
to the noise level faster than the y power. The noise filter seems 
to be correct in that it removes the noise without visibly chang- 
ing the resolution in the smallest structures. This is also the case 
for the noise filters constructed for the other tilt angles. 

For a proper Wiener-type noise filter, one needs the sig- 
nal power and the noise power. The original, raw-data noise 



(a) 6?4 



(b) 4°8 



(c) 4?2 



Fig. 15: The low-pass noise filter for the largest angles. The 
perimeter of the circle represents the diffraction limit. 



at the highest spatial frequencies is already filtered out by the 
MOMFBD restoratior^ Instead, there is a noise that comes from 



^ In the raw data, assuming additive Gaussian noise, the MOMFBD 
program measures the noise level outside the diffraction limit where we 
know there is no signal. 
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Table 1: Observed RMS contrasts. 



Deconvolution 








NB 








WB 


r.2 


2°5 


3?0 


3?7 


4?2 


4°8 


6°4 


^ atm 

S atmS tilt(^) 


25.7 
25.8 


16.5 
16.6 


15.4 
15.6 


13.1 
13.6 


12.0 
12.4 


12.2 
12.7 


12.9 
13.6 


13.3 



Notes. RMS contrasts of real NB data in percent of the mean intensity. The first row has contrasts for images corrected for atmospheric turbulence 
by the MOMFBD program. The second row contrasts are for the same images after compensation also for Smt. WB contrasts are given for 
comparison with the 6° 4 data. 



the mosaicking of subfields done as part of the MOMFBD pro- 
cessing. The power of this noise seems to vary with spatial fre- 
quency in a way similar to the signal but extend to higher spatial 
frequencies. The filters in Fig. 15 are based on thresholding the 
power spectrum of the uncompensated NB image at levels found 
by trial and error. The filters are then made from the resulting bi- 
nary images by closing holes and removing isolated pixels, and 
finally smoothing with a Gaussian kernel. 

For images to be used for making diflTerential quantities (like 
magnetograms or Doppler maps), it is best to use the same low- 
pass noise filter. Otherwise, there will be spatial frequencies with 
information missing in some images but not in others. A simple 
way of constructing such a filter is to form, in every Fourier do- 
main pixel, the minimum of the filters for the relevant images. 
Such a combination of the filters shown in Fig. 15 would be very 
similar to the 6? 4 filter. 

The contrasts in Table [T] are much lower than those in Fig.|9] 
because of various sources of stray light, see discussion by 
|S char mer et al.| ( [20TQ| . The contrasts show the expected increase 
in contrast with Stnt compensation, more for large 6 than for 
small. 



5.3. Image shift and MOMFBD processing 

A very useful property of the Multi-Object part of MOMFBD 
image restoration is that the restored images of the diflTerent ob- 
j ects are delivered by th e program, aligned to subpixel precision 
( [van Noort et al.|20^6] ). In this context, different "objects" refer 
to the same patch of the Sun, but imaged in different wavelengths 
and/or different polarization states. This ability facilitates the 
calculation of physical quantities based on combinations of these 
objects with a minimum of artifacts from misalignment. When 
the images are co-spatial and co-temporal, they are subject to 
the same wavefront aberrations and are therefore blurred and 
shifted by the same PSFs (or PSFs that differ in easily mod- 
eled ways). The alignment is accomplished by a pre-processing 
step, where local shifts are measured in images of an array of 
pinholes, which provides enough information for the MOMFBD 
program to compensate for fixed differences in alignment and, 
to some extent, in rotation or image scale. 

Compensation for 5'tiit interferes with this procedure because 
the corresponding PSF is off-center and therefore shifts the im- 
ages. If the MOMFBD-restored images are shifted in a post- 
processing step, it is necessary to make sure the pinhole array 
images are shifted the same way before the local- shift measure- 
ments are done. So the pinhole array images should be 5'tiit com- 
pensated before the pre-processing step. 

However, the difference in shifts for large and small tilt an- 
gles can be several pixels, see Fig. |7] Applying the compensa- 
tion to the pinhole images means the subfield grid used for the 
corresponding objects will be off by the same number of pix- 
els. This is then corrected for in the post-processing step, where 



the restored images are 5'tiit compensated, but it is likely that 
the alignment works better if the subfield grids are aligned dur- 
ing the MOMFBD process. This can be accomplished if the 5'tiit 
PSFs are centered before they are used. We believe it is enough 
and possibly safest to do this to nearest pixel and avoid subpixel 
operations for this step. 

We therefore recommend the following steps, which we have 
tested with our data set: 

1. Calculate 5'tiit(^) for all 6 and center the corresponding PSFs 
to nearest pixel precision. 

2. Deconvolve the pinhole array images with these centered 
PSFs. 

3. Perform the standard pre-processing step using the decon- 
volved pinhole array images. 

4. Do normal MOMFBD processing. 

5. Deconvolve the restored images with the centered PSFs. 



6. Conclusions 

We have investigated a NB double-cavity filter, used at the SST 
for scanning through the blue wing of the Can H spectral line. 

Measurements with the TRIPPEL spectrograph show that 
the filter has a passband that is narrower by 10% than indicated 
by profiles obtained from the manufacturer. We note that the fil- 
ter was 8 years old at the time of these measurements but we 
know too little about coating age effects to speculate on why this 
can happen. These measurements also indicate that a descrip- 
tion based on wavelength shifts of single-cavity filter profiles, 
together with a model for variations in the incident angle over 
the pupil, work very well for predicting the position and shape 
of the broadened passband. 

We calculate PSFs and OTFs corresponding to a diffraction 
limited telescope and tilt angles in the range 0?0-7?5 and find 
that the Strehl ratio decreases to less than 0.6 for the largest an- 
gle. 

We find that the tilt effects couple very weakly with wave- 
front corrections and that the OTFs are separable in the usual 
OTF for phase aberrations (including atmospheric turbulence ef- 
fects) and 5'tiit, a correction for the filter tilt effects. This means 
it is sufficient to correct images after restoration for atmospheric 
effects, the raw data do not have to be corrected before im- 
age restoration. However, the PSFs are in general off-center and 
asymmetric, which has consequences for MOMFBD alignment. 
We have designed a procedure for 5'tiit compensation in a way 
that maintains the alignment properties of MOMFBD image 
restoration. 

By experiments with synthetic data, we find that deconvo- 
lution with the proper OTF compensates very well for the de- 
grading effects that are really affecting the image formation at 
each wavelength in the passband separately. We demonstrate the 
effect of the phase of the transmission profile and show that not 
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Fig. A.l: Filter tilt geometry, (a) Projected on the pupil plane, 
(b) Top view, projection on the horizontal plane. The pupil with 
diameter D = lis indicated by the circle in a) and by the two tick 
marks on the x axis in b). The focal length is /. In (x,y, z) space 
(z is the optical axis), a ray leaves the pupil plane at (Xp,3;p,0) 
and hits the detector in the focal plane at (xf,yf,F). The filter 
surface normal from (xf,yf, F) intersects with the pupil pl ane a t 
an x-coordinate increased by an amount Ax given by Eq. ( |A.5| ). 
The angle ^ is in a horizontal plane but a is in general not. 

allowing for this phase in the right way can significantly increase 
the RMS error in the compensated images. 

Our theoretical model for pupil apodisation eflTects from 
varying incident angles on NB interference filters does not ade- 
quately predict the eflTects on image quality. This is evident from 
power spectra of pinhole images as well as solar data. This can 
be an error in our derivations or coding, but we note that it would 
have to be a mistake that aflfects 5'tiit without causing errors in the 
profile broadening. 

Instead, we find that useful 5'tiit can be measured by use of 
pinhole calibration data. The ones we measure are visibly af- 
fected by noise but this can easily be improved by use of data 
with better SNR. Note that the better the solar data, the better 
pinhole images are needed. 

For the specific setup with a tilt-tunable Ca ii H filter at the 
SST, the degrading effect on the PSF is negligible when the filter 
is not tilted, but becomes more significant as the tilt angle is 
increased. 

A filter wheel with fixed passband filters or an FPI are better 
ways of sampling diflferent positions in a spectral line. However, 
tuning interference filters by tilting them from normal incidence 
is useful because it is a cheaper solution. It may be preferable 
if one can live with the spectral broadening and compensate for 
the degraded resolution. 

Acknowledgements. Goran Scharmer initiated this research severa l years ago, 
early results were presented already by |Lofdahl & Scharmer| ( |2004) , and he has 



also provided many useful comments after the project was recently revived. Luc 
Rouppe van der Voort took part in the observational testing of the filter with 
TRIPPEL. Mats Carlsson is thanked for sharing the 3D MHD simulation snap- 
shot used for making synthetic data. The Swedish 1-m Solar Telescope is oper- 
ated on the island of La Palma by the Institute for Solar Physics of the Royal 
Swedish Academy of Sciences in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 



References 

Beckers, J. M. 1998, A&AS, 120, 191 

Brault, J. W. & Neckel, H. 1987, Spectral Atlas of Solar Absolute Disk- 
Averaged and Disk-Center Intensity from 3290 TO 12510 A, available from 
Hamburg Observatory Anonymous FTP Site since 1999. See announcement 
in Sol. Phys. 184, 421 

Carlsson, M. 1986, A computer program for solving multi-level non-LTE radia- 
tive transfer problems in moving or static atmospheres, Uppsala Astronomical 
Observatory Reports 33 

Kiselman, D. 2008, Physica Scripta, T133, 014016 

Klein, M. V. & Furtak, T. E. 1986, Optics, 2nd edn. (New York: John Wiley & 
sons. Inc.) 

Kupka, R, Piskunov, N., Ryabchikova, T. A., Stempels, H. C, & Weiss, W. W. 

1999, A&AS, 138, 119 
Kupka, F. G, Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C, & Weiss, 

W. W. 2000, Baltic Astronomy, 9, 590 
Lofdahl, M. G. 2002, in Proc. SPIE, Vol. 4792, Image Reconstruction from 

Incomplete Data II, ed. P. J. Bones, M. A. Fiddy, & R. P. Millane, 146-155 
Lofdahl, M. G. & Scharmer, G. B. 2004, in Solar Image Processing Workshop II 

(abstracts), ed. P. Gallagher, Annapolis, Maryland, USA 
Noll, R. J. 1976, J. Opt. Soc. Am., 66, 207 

Piskunov, N. E., Kupka, R, Ryabchikova, T. A., Weiss, W. W., & JefFery, C. S. 

1995, A&AS, 112, 525 
Potter, J. R. 2004, private communication by email, response to web 

query to Barr Associates, Inc. in December 2004 via htt p : //www . | 

barrassociates . com/contact . php 
Press, W. H., Flannery, B. R, Teukolsky, S. A., & Vetterling, W. T. 1986, 

Numerical Recipes, The Art of Scientific Computing (Cambridge University 

Press) 

Rouppe van der Voort, L. 2002, PhD thesis, Stockholm University 
Scharmer, G. B. 2006, A&A, 447, 1111 

Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Pettersson, B. 

2003a, in Proc. SPIE, Vol. 4853, Innovative Telescopes and Instrumentation 

for Solar Astrophysics, ed. S. Keil & S. Avakyan, 341-350 
Scharmer, G. B., Dettori, R, Lofdahl, M. G, & Shand, M. 2003b, in Proc. SPIE, 

Vol. 4853, Innovative Telescopes and Instrumentation for Solar Astrophysics, 

ed. S. Keil & S. Avakyan, 370-380 
Scharmer, G. B., Lofdahl, M. G, van Werkhoven, T. I. M., & de la 

Cruz Rodriguez, J. 2010, A&A, 521, A68 
Smith, W. J. 1990, Modern Optical Engineering: The Design of Optical Systems, 

2nd edn. (New York: McGraw-Hill) 
Stein, R. R & Nordlund, A. 1998, ApJ, 499, 914 

van Noort, M., Rouppe van der Voort, L., & Lofdahl, M. G. 2005, Sol. Phys., 
228, 191 

von der Luhe, O. & Kentischer, T. J. 2000, A&AS, 146, 499 



Appendix A: The incident angle 

This appendix gives a derivation of an expression for the incident 
angle a on a filter in a converging beam as a function of position 
in the pupil. We use a coordinate system with the unit length 
equal to the pupil diameter, D. The distance between the pupil 
and focal planes is then equal to the F-ratio, F = f/D, where / 
is the focal length. 

A ray emanates from the pupil plane at the point (x^^y^, 0) 
and intersects the detector in the focal plane at (xf,yf,F). The 
filter can be anywhere between the pupil and the detector and it 
is tilted by an angle 6 around an axis parallel to the y axis. The 



geometry of the pupil-filter system is shown in Fig. A.l 



The Law of Cosines gives a relation between the incident 
angle a and the lengths of the sides in the abc triangle, 



2bc cos a. 



(A.1) 
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where a, /?, and c are the distances between the triangle corners, 
= (Xf + Ax - x^f + (jf - 3;p)^ (A.2) 

b^ = {x^-Xff + {y^-yif + F\ (A.3) 
= + (A.4) 

The length Ax is given by the filter tilt angle 6 as 

Ax = F tan a (A.5) 



Mats G. Lofdahl et al.: A tilted interference filter in a converging beam 

For 



cos Of : 



<7il2, Eqs. ( |A.1| )-( [A3] ) and trivial algebra yield 
F cos ^ + (Xp - Xf ) sin 6 



(A.6) 
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